SIMU=0;CD_REAL_500_0=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_0[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_0=c(CD_REAL_500_0,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_0,'CD_REAL_500_0.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_1=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_1[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_1=c(CD_REAL_500_1,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_1,'CD_REAL_500_1.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_2=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_2[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_2=c(CD_REAL_500_2,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_2,'CD_REAL_500_2.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_3=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_3[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_3=c(CD_REAL_500_3,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_3,'CD_REAL_500_3.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_4=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_4[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_4=c(CD_REAL_500_4,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_4,'CD_REAL_500_4.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_5=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_5[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_5=c(CD_REAL_500_5,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_5,'CD_REAL_500_5.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_6=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_6[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_6=c(CD_REAL_500_6,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_6,'CD_REAL_500_6.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_7=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_7[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_7=c(CD_REAL_500_7,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_7,'CD_REAL_500_7.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_8=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_8[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_8=c(CD_REAL_500_8,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_8,'CD_REAL_500_8.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_9=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_9[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_9=c(CD_REAL_500_9,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_9,'CD_REAL_500_9.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_10=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_10[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_10=c(CD_REAL_500_10,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_10,'CD_REAL_500_10.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_15=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_15[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_15=c(CD_REAL_500_15,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_15,'CD_REAL_500_15.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_20=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_20[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_20=c(CD_REAL_500_20,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_20,'CD_REAL_500_20.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_25=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_25[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_25=c(CD_REAL_500_25,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_25,'CD_REAL_500_25.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_30=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_30[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_30=c(CD_REAL_500_30,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_30,'CD_REAL_500_30.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_35=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_35[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_35=c(CD_REAL_500_35,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_35,'CD_REAL_500_35.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_40=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_40[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_40=c(CD_REAL_500_40,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_40,'CD_REAL_500_40.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_45=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_45[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_45=c(CD_REAL_500_45,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_45,'CD_REAL_500_45.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_50=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_50[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_50=c(CD_REAL_500_50,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_50,'CD_REAL_500_50.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_55=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_55[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_55=c(CD_REAL_500_55,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_55,'CD_REAL_500_55.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_60=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_60[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_60=c(CD_REAL_500_60,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_60,'CD_REAL_500_60.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_65=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_65[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_65=c(CD_REAL_500_65,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_65,'CD_REAL_500_65.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_70=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_70[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_70=c(CD_REAL_500_70,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_70,'CD_REAL_500_70.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_75=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_75[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_75=c(CD_REAL_500_75,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_75,'CD_REAL_500_75.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_80=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_80[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_80=c(CD_REAL_500_80,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_80,'CD_REAL_500_80.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_85=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_85[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_85=c(CD_REAL_500_85,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_85,'CD_REAL_500_85.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_90=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_90[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_90=c(CD_REAL_500_90,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_90,'CD_REAL_500_90.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_95=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_95[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_95=c(CD_REAL_500_95,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_95,'CD_REAL_500_95.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_100=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_100[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_100=c(CD_REAL_500_100,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_100,'CD_REAL_500_100.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_110=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_110[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_110=c(CD_REAL_500_110,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_110,'CD_REAL_500_110.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_120=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_120[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_120=c(CD_REAL_500_120,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_120,'CD_REAL_500_120.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_130=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_130[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_130=c(CD_REAL_500_130,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_130,'CD_REAL_500_130.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_140=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_140[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_140=c(CD_REAL_500_140,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_140,'CD_REAL_500_140.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_150=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_150[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_150=c(CD_REAL_500_150,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_150,'CD_REAL_500_150.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_160=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_160[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_160=c(CD_REAL_500_160,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_160,'CD_REAL_500_160.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_170=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_170[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_170=c(CD_REAL_500_170,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_170,'CD_REAL_500_170.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_180=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_180[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_180=c(CD_REAL_500_180,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_180,'CD_REAL_500_180.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_190=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_190[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_190=c(CD_REAL_500_190,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_190,'CD_REAL_500_190.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_200=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_200[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_200=c(CD_REAL_500_200,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_200,'CD_REAL_500_200.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_210=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_210[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_210=c(CD_REAL_500_210,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_210,'CD_REAL_500_210.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_220=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_220[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_220=c(CD_REAL_500_220,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_220,'CD_REAL_500_220.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_230=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_230[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_230=c(CD_REAL_500_230,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_230,'CD_REAL_500_230.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_240=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_240[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_240=c(CD_REAL_500_240,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_240,'CD_REAL_500_240.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_250=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_250[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_250=c(CD_REAL_500_250,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_250,'CD_REAL_500_250.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_260=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_260[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_260=c(CD_REAL_500_260,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_260,'CD_REAL_500_260.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_270=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_270[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_270=c(CD_REAL_500_270,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_270,'CD_REAL_500_270.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_280=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_280[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_280=c(CD_REAL_500_280,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_280,'CD_REAL_500_280.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_290=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_290[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_290=c(CD_REAL_500_290,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_290,'CD_REAL_500_290.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_300=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_300[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_300=c(CD_REAL_500_300,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_300,'CD_REAL_500_300.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_310=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_310[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_310=c(CD_REAL_500_310,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_310,'CD_REAL_500_310.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_320=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_320[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_320=c(CD_REAL_500_320,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_320,'CD_REAL_500_320.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_330=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_330[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_330=c(CD_REAL_500_330,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_330,'CD_REAL_500_330.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_340=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_340[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_340=c(CD_REAL_500_340,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_340,'CD_REAL_500_340.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_350=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_350[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_350=c(CD_REAL_500_350,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_350,'CD_REAL_500_350.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_360=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_360[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_360=c(CD_REAL_500_360,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_360,'CD_REAL_500_360.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_370=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_370[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_370=c(CD_REAL_500_370,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_370,'CD_REAL_500_370.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_380=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_380[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_380=c(CD_REAL_500_380,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_380,'CD_REAL_500_380.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_390=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_390[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_390=c(CD_REAL_500_390,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_390,'CD_REAL_500_390.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_400=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_400[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_400=c(CD_REAL_500_400,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_400,'CD_REAL_500_400.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_410=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_410[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_410=c(CD_REAL_500_410,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_410,'CD_REAL_500_410.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_420=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_420[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_420=c(CD_REAL_500_420,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_420,'CD_REAL_500_420.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_430=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_430[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_430=c(CD_REAL_500_430,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_430,'CD_REAL_500_430.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_440=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_440[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_440=c(CD_REAL_500_440,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_440,'CD_REAL_500_440.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_450=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_450[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_450=c(CD_REAL_500_450,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_450,'CD_REAL_500_450.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_460=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_460[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_460=c(CD_REAL_500_460,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_460,'CD_REAL_500_460.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_470=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_470[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_470=c(CD_REAL_500_470,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_470,'CD_REAL_500_470.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_480=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_480[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_480=c(CD_REAL_500_480,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_480,'CD_REAL_500_480.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_490=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_490[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_490=c(CD_REAL_500_490,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_490,'CD_REAL_500_490.txt',row.names=FALSE,col.names=FALSE)

SIMU=0;CD_REAL_500_500=numeric()
for(w in 1:10000)
{
y1=PROFILES_G1_500[,w];y2=PROFILES_G2_500[,w]
DIF=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
dif=(x-1)*stap-(y-1)*stap
DIF[x,y]=dif
}}
DIF=sign(DIF)
DEN=matrix(0,length(y1),length(y2))
for(x in 1:length(y1))
{
for(y in 1:length(y2))
{
den=y1[x]*y2[y]*stap^2
DEN[x,y]=den
}}
CD_real=sum(DIF*DEN)/sum(DEN)
CD_REAL_500_500=c(CD_REAL_500_500,CD_real)
SIMU=SIMU+1;if(SIMU%%1000==0)print(SIMU)
}
write.table(CD_REAL_500_500,'CD_REAL_500_500.txt',row.names=FALSE,col.names=FALSE)

